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Abstract 

We have studied galaxy two-point angular correlations as a function of color using 
4-m plate photometry in two independent fields. Each field consists of over 2900 galaxies 
with magnitudes 20 < Bj < 23.5 in an area of approximately 750 arcmin 2 . We find 
that the autocorrelation amplitude of the bluest 15% of galaxies is surprisingly strong, 
with a relative increase in clustering amplitude of a factor of 6 over that of the complete 
data set, while exhibiting a power law slope consistent with the canonical value of —0.8. 
These very blue galaxies are also found to be weakly correlated with galaxies of median 
color and marginally ant i- correlated with the reddest subset. These correlation properties 
are incompatible with existing simple models of the galaxy distribution; they suggest that 
a significant fraction, more than 50%, of these very blue galaxies are a faint population 
which lie at nearby redshifts z < 0.3. 
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1. Introduction 



The nature of faint field galaxies and their high surface density at faint magnitudes 
remain a mystery despite extensive new observations over the last 15 years (see Koo and 
Kron 1992, Lilly 1993, and Koo 1996 for recent reviews). Besides fainter number counts 
and colors, and more recently the morphology of faint galaxies from Hubble Space Tele- 
scope observations (Griffiths et al. 1994), the key new observations have been redshifts to 
limits as faint as Bj = 24 for 12 galaxies (Cowie, Songaila, & Hu 1991), and 73 galaxies 
(Glazebrook et al. 1995) and to / = 22 for almost 600 galaxies (Lilly et al. 1995). The 
detailed conclusions of these deep redshift surveys remain at odds, though all agree that 
substantial evolution has occurred in the luminosity function of field galaxies (Lilly et al. 
1995, Colless 1995). Whether such changes are caused by a new population, substantial 
mergers, or enhanced star-formation among dwarfs has yet to be resolved. Moreover, both 
of the larger surveys remain incomplete at the 20% to 27% level, with the bulk of the 
unidentified galaxies being bluer than average. This level of incompleteness is high enough 
to make a significant difference in understanding the nature of faint galaxies, especially if 
the bulk of the unidentified population is either at high or low redshifts. 

This paper places additional constraints on the nature of the bluest galaxies by using 
another observable property of faint galaxies, namely, their angular clustering behavior. 
Although the approach of using angular correlation functions to study galaxy evolution is 
not new (Koo & Szalay 1984) and has been applied to several samples of faint galaxies, 
at different wavelengths, or for gross cuts in color (Efstathiou et al. 1991, Bernstein et al. 
1994, Infante & Pritchet 1995, Neuschaefer & Windhorst 1995), none of the analyses so far 
have explored the variation of clustering among galaxies of extreme colors. Here we use 
the angular correlation function as a tool to measure, or at least constrain, the redshift 
distribution of such galaxies. 

The two-point angular correlation function w(9) is well known and defined by the 
joint probability 5P of finding galaxies in both of the elements of solid angle 5Q± and 5Q 2 
placed at an angular separation 9, 

5P = N 2 [l + w(e)]5tt 1 5tt 2 (1) 

where N is the mean surface density of galaxies. The observations are usually close to a 
power law, 

w(9) = A w 6 x -\ (2) 

with amplitude A w and 7 ~ 1.8 (see Peebles 1980). As an estimator for this function we 
will use that described in Landy & Szalay (1993) for its improved error properties. 

The scaling of angular correlations as a function of limiting magnitude is well under- 
stood: as the limiting magnitude of a survey increases, there are more projected pairs that 
weakens the correlations (Peebles 1980), and the amplitude is proportional to an integral 
over the second moment of dn/dz, the redshift distribution of the galaxies, 



A w oc I dz c(z) [ 
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where c(z) depends on evolution, curvature, and 7 (see e.g. Koo & Szalay 1984). Less is 
known about changes due to color selection. Differences in the measured clustering as a 
function of angle w(0), for subsamples such as those chosen by color, can occur for two 
reasons: each subsample may be dominated by a different population with different cluster- 
ing properties, and/or different redshift distributions. This degeneracy between intrinsic 
clustering and redshift distributions can be broken by considering other information such 
as the overall shape of the correlation function, cross-correlations with known populations, 
and number counts. Correlation functions are normally not thought of as a probe of red- 
shift distributions, but as pointed out earlier (Koo and Szalay 1984), angular correlation 
functions provide powerful constraints on this distribution (and hence evolution and cos- 
mology), which is particularly valuable in the absence of spectroscopic redshifts for very 
faint samples, for example, to / > 26. Thus the method serves as an independent check or 
constraint on models that attempt to describe the more easily measured colors and counts 
of faint galaxies. 

Since early- type (E/SO/Sa) galaxies have stronger than average intrinsic clustering 
(Giovanelli, Haynes, & Chincarini 1986), the a priori expectation is for red galaxies to 
have correlation amplitudes higher than that of an average field sample. In contrast, 
late-type, blue galaxies today appear to be more weakly clustered than an average field 
sample (Nicotra, Abadi, & Lambas 1993). Some groups have proposed that such galaxies 
were brighter in the past because of enhanced star formation (Broadhurst, Ellis, & Shanks 
1988). If so, their redshift distributions will be more extended than without this luminosity 
increase. In this case, one might expect a further reduction in clustering amplitude for 
blue galaxies observed to fainter limiting magnitudes. It has also been recently reported 
that low surface brightness galaxies are more weakly clustered than the general population 
(Neuschaefer & Windhorst 1995), also expected in a CDM dominated universe (Efstathiou 
1995). 

To compare these expectations with observations, we have calculated the correlation 
functions of subsets of galaxies chosen as a function of color: U — Bj, U — Rp, and Bj — R F . 
We start with extreme subsets of either red or blue objects and then increase the range 
to include objects of less extreme color. Also calculated is the correlation function for 
samples selected in bands of color as well as the cross correlation between samples selected 
in different bands. We report results on U — Rp, as this color gives the largest wavelength 
range to discriminate color-dependent effects. In addition, this is the combination for which 
the galaxy colors are least explicitly dependent on their redshift, and thus most correlated 
with intrinsic color (Koo 1986). The correlation amplitude dependence on U — Bj and 
Bj — Rp is similar to U — Rf, although the enhancements were weaker. 

2. Data 

The photometry data are from two fields of approximately 750 arcmin 2 each and 
based on photographic plates taken at the prime focus of the Mayall 4m telescope at Kitt 
Peak National Observatory. One field SA57 (1306+2939) is at the North Galactic Pole; 
the other field SA68 (0015+1537) has galactic coordinates of I = —111 and b = —46 and 
is near Perseus-Pisces. Details of the observations and data reduction techniques have 
been published elsewhere (Kron 1980, Koo 1986). Here it is sufficient to note that the 
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plates were digitized using the Berkeley PDS microdensitometer; that areas in the fields 
around bright stars as well as bad pixel lines in the data reduction have been excised from 
the fields; that both fields are considered complete to Bj = 23.5, and finally that two 
plates were used in each bandpass with only multiple detections included in the catalogs 
to decrease the incidence of artifacts. We have selected galaxies in the magnitude range 
20 < Bj < 23.5, with 2922 galaxies in SA57 and 2999 galaxies in SA68. 

The correlation properties of the same data have previously been analyzed, but only 
as a function of Bj magnitude (Koo & Szalay 1984). This earlier work is in agreement with 
that of other researchers (Stevenson et al. 1985; Pritchet & Infante 1992; Jones, Shanks, 
& Fong 1988; Bernstein et al. 1994; Brainerd, Smail, & Mould 1995 ), who found that 
the scaled two-point angular correlation function amplitude is consistent with that found 
locally, except for hints of mild evolution in luminosity and/or clustering. 

Given the unexpected nature of our results (see below) we have performed numerous 
checks to determine the validity of our measurements. It should, however, be emphasized 
that these data are derived from multiple plates and that the results are qualitatively 
identical for the three independent color combinations for both red and blue galaxies 
in two separate fields. The similarity of results decreases the probability that artifacts, 
such as gradients across the plates in different passbands, are responsible for the signal. 
As a further check against artifacts caused by gross variations in color over the plates, 
the correlation analyses were rerun with the median U — Rp color equilibrated over the 
plates by use of a 400"x400" square smoothing window. Although this adjustment, which 
amounted to several tenths of magnitude in some areas, weakened the signals, all general 
trends were still clearly evident. Additional tests were run on the data to check whether 
either population was associated with bright objects such as galaxies or stars in the fields. 
These proved negative. Cross correlations between stars and galaxies in the fields also 
indicated no obvious problems with the catalogs as did analysis using a subset of stars 
identified through proper motion analysis. 

3. Results 

3.1 Measured Angular Correlation Amplitudes 

The correlation amplitudes for SA57 and SA68, as a function of U — Rf color, are 
shown in Figure 1, together with plots showing the percentage of data included in each 
subset. These results show the angular correlation amplitude of the reddest and bluest 
subsets of galaxies, respectively. The horizontal line is the correlation amplitude of the 
whole sample. In order to present a single number, an amplitude, the correlation function 
was integrated between 9 and 144 arcsecs to reduce noise and for compatibility with the 
earlier work (Koo & Szalay 1984). Assuming a nominal redshift of 0.3, H Q = 50 km 
sec _1 Mpc _1 , and q Q ~ these limits span 70 kpc to 1.1 Mpc. The amplitude of w(6), 
A w , is reported using the value of —0.8 for the power law slope at one degree in units 
of 10~ 4 . For the average of the two fields we obtained A w (all) = 18.0 ± 1.8. Statistical 
(Poisson) error bars are shown for all results, although these certainly underestimate the 
true variance in the signal, as will be discussed below. 

Both fields show a marked enhancement in correlation amplitude, with increases 
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greater than a factor of ten for the extreme subsets in both red and blue. As a check 
against the possibility that only the most extreme galaxies contribute to the rise in ampli- 
tude (because of artifacts in the data), we calculated w{9) in bands of color for both the 
red and blue galaxies. The results are shown in Figure 2. It is evident that the increase in 
angular correlation function exists over more than one magnitude in color, which indicates 
that the enhancement is genuine. 

Figure 3 shows w{6) for 'red' galaxies U—Rf > 1.5 and 'blue' galaxies U—Rf < —0.25 
together with the entire data set for both fields in five bins of equal logarithmic width 
between 18" and 289" . These results are compared against a fiducial no-evolution model 
assuming a power law index of -0.8, q o =0- 01, and H o =50 km/sec. Both the red and blue 
subsets in each field show qualitative agreement with a power law slope of -0.8, together 
with enhanced clustering over this entire range. 

Other recent correlation function studies (Efstathiou et al. 1991; Brainerd et al. 1995; 
Neuschaefer et al. 1995) that used CCD detectors have found quite weak angular correlation 
functions for objects with an isophotal magnitude 24 < Bj < 26. When the Efstathiou 
et al. (1991) sample was divided into red and blue subsets, the blue subset was found 
to have a marginally (< 2a) stronger correlation amplitude. One of their fields, SA68, 
overlaps with one of ours. When we apply a similar color split to our data as adopted by 
these authors, we obtain similarly low amplitudes (as seen on Figure 1). Only with more 
extreme color cuts do the enhancements in correlation amplitude become evident. 

Our error bars correspond to the statistical noise in the LS estimator (Landy & Szalay 
1993), i.e. we assume that the distribution of the total number of galaxies in the field 
is Poisson, thus the variance in the correlation is Poisson in the bin counts. The LS 
estimator, conditional on the number of galaxies observed, is optimal in this case, and still 
minimizes the variance in all the other cases. If, however, there is additional large scale 
clustering (as shown by the presence of a cluster), the irreducible higher order moments of 
the galaxy counts on scales of the whole field will result in additional variance. An a priori 
determination of this variance, on the other hand, requires knowledge of these higher order 
moments at these magnitudes, colors, and angular scales. As these are not available, we 
only plot the Poisson errors, but expect the full error to be about a factor of two higher. 

Table 1 

Auto and Cross- Correlations by Bands of Color 

Blue (20%) Mid (60%) Red (20%) 

Blue 68± 20 9± 7 -18± 10 

Mid 14± 1 20± 9 

Red 130± 30 



Table 1 reports the auto- and cross-correlations of the top 20% blue and red galaxies 
along with the mid 60 %. The error bounds in the table are derived from the actual scatter 
between the two fields, and appear to be twice the statistical, Poisson errors. These results 
show that the auto-correlations of the intermediate colors are weaker than for the whole 
data set, and the strong correlations in the two extreme colors boost the overall clustering 
amplitude to its observed value 18.0 ± 1.8. In addition, the cross correlation function 



5 



between the red and blue populations is marginally anti-correlated. 

3.2 Comparison with Models 

We have compared our results to predictions of models (Gronwall & Koo 1995), de- 
signed to fit the observed counts, colors, and redshift distributions, but that rely on having 
the luminosity function of field galaxies be a free "parameter", i.e., one that is not con- 
strained by local observations. The models include some mild evolution in the stellar 
populations and a small amount of intrinsic absorption/reddening due to dust. The use 
of models was not meant to be part of an exhaustive analysis of various evolutionary and 
cosmological scenarios, but rather to provide vital baseline estimates of the redshift distri- 
butions needed to interpret the angular correlations. These models are improved estimates 
over those of prior analyses of angular correlation functions (Efstathiou et al 1991, Koo 
& Szalay 1984), since they are also found to fit a wide range of other faint galaxy counts, 
colors, and redshifts. 

To compare the actual observations, simulated galaxy catalogs were generated that 
account for photometric errors and incompleteness that mimic those expected to exist 
in our data. We then applied identical magnitude (20 < Bj < 23.5) and color cuts 
(U — Rf > 1-5 and U — Rf < —0.25) to the simulated catalogs and derived the redshift 
distributions. These were used to calculate the angular correlation amplitude, assuming a 
standard no-evolution, -1.8 slope correlation function. The models reproduce the overall 
correlation amplitude accurately (18.2), and can be reconciled with the amplifications seen 
in the 'red' subsets, but fail to explain the angular clustering of the extremely 'blue' subset, 

A w (b\ue)/A w (al\) =0.44 

A w (red)/A w (all) =1.67 (4) 
A w (cross) /A w (all) =1.02 

The discrepancy is so large that no minor modification of the models can give a satisfactory 
agreement. The redshift distribution for the 'blue' galaxies was more extended than for 
the whole sample (because of evolution), resulting in a larger number of projected pairs, 
thus in a lower amplitude for the predicted angular correlation amplitude. Thus the main 
question is what constraints these observations generate on the redshift distribution of the 
bluest galaxies. 

3.3 Comparison to Redshift Surveys 

Deep redshift surveys are just reaching our magnitude limits (Colless et al. 1991, 1993; 
Glazebrook et al. 1995; Lilly et al. 1995). The redshift survey most similar in limiting 
magnitude to our data is the Glazebrook et al. (1995) survey to B < 24. This sample 
has 73 measured redshifts. Taking the 21 bluest galaxies of the B < 23.5 subset, only 10 
have redshifts; thus the incompleteness is over 50%, whereas the whole survey was ~ 70% 
complete. These 10 galaxies were rather evenly distributed over the redshift range of 0.08 
to 1.1; thus this blue subset of the spectroscopic sample is too small and incomplete to draw 
reliable conclusions about the detailed redshift distribution to be drawn when such narrow 
color divisions are applied. The Lilly et al. (1995) sample to / < 22 has a larger number 
of redshifts (600), but these are unpublished. We discuss their result in the Summary. 
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4. Discussion and Analysis 
4-1 The Red Galaxies 

In some respects, the large enhancement in the correlation amplitude of the reddest 
galaxies may seem as surprising as that for the blue subset. However, this result can be 
explained in terms of known properties of red galaxies quite naturally. An enhancement in 
the correlation amplitude for the red galaxies may be expected since these are primarily 
early type (E/SO) galaxies (morphology-density relation). Further inspection of the red 
subsets found in these fields reveals that field SA57 contains a known cluster at z = 0.24 
and SA68 contains a supercluster at a redshift of z = 0.54, both visible in maps of the fields 
(Figs. 10 and 11 in Koo 1986). From the models we have estimated that, for galaxies redder 
than U — Rf > 1-5 (with a median redshift of 0.34) the correlation amplitude is enhanced 
by a factor 1.67 from the modified redshift distribution alone. Since that population is 
hardly evolving, these model predictions should be quite robust. On the other hand, a 4 
times stronger intrinsic clustering has actually been measured for E/SO galaxies (Nicotra 
et al. 1993), resulting in a total correlation amplification of 6.4. Therefore, the measured 
overall enhancement of the red galaxy signal can be well understood in these terms, and 
illustrates the sensitivity of using angular correlations as a function of color to study both 
galaxy evolution and large scale structure. 

4-2 Blue Galaxies 

Analyzing the enhancement in the correlation amplitude of the bluest galaxies is more 
problematic than that of the red subset since these have not been identified securely with 
a well-studied local population. To reach quantitative conclusions about their redshift 
distribution, we will exploit two basic observational constraints on blue galaxies which 
generate opposing trends: (1) the amplification of the very blue galaxy angular correlations 
is greater than 6, and (2) the population fraction is greater than 15%. These constraints 
provide strong limits on the possible redshift distributions. 

Several hypotheses have been proposed to explain the nature and distribution of faint 
blue galaxies. One hypothesis is that faint blue galaxies are the precursors of normal 
galaxies but undergoing enhanced star formation at an earlier epoch, z > 1. Another is 
that these galaxies are a merging population, perhaps at intermediate to high redshifts, 
since merging galaxies are expected to be fairly bright and blue because of enhanced star 
formation (Broadhurst, Ellis, & Glazebrook 1992). Other hypotheses include a dominant 
population of relatively nearby, low surface brightness galaxies that are not represented in 
existing estimates of the local luminosity function (McGaugh 1994) and the existence of a 
local bursting dwarf population with strong luminosity evolution that has presently faded 
from view (Broadhurst et al. 1988; Babul & Rees 1992; Eales 1993; Lilly 1993). Below, we 
discuss these possibilities to assess whether any of these scenarios is consistent with our 
correlation function measurements. 

An enhancement in autocorrelation amplitude for the blue subsample can be explained 
either by these galaxies being intrinsically more strongly correlated than commonly as- 
sumed, or their redshift distribution being very different, or a combination of both. Given 
that the intrinsic spatial correlation of late-type galaxies has been found to be 4 times 
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weaker than that of early type galaxies (Nicotra et al 1993), arguments that postulate 
increased correlation amplitudes appear to be untenable. To explain the observed enhance- 
ment in w(9) solely by this effect would require the opposite trend, namely, a 12 times 
stronger intrinsic clustering for the very blue subset if they had the redshift distribution 
as predicted by the model, or 6 times stronger if they had the redshift distribution of 
the whole sample. This leaves the possibility that the redshift distribution of the bluest 
galaxies is quite different from the predictions. Since the angular correlation amplitude is 
proportional to the second moment of dn/dz, the most obvious implication of the enhanced 
clustering amplitude is that the redshift distribution is narrower than the models suggest. 

In considering the possibility that the bluest galaxies are all at high redshift, let us 
assume that our blue galaxies are the distant counterparts of normal galaxies today that 
have undergone enhanced star formation in the past. Having measured the power-law 
slope of the correlation function of the blue galaxies to be similar to that of galaxies today, 
and assuming that these counterparts predominantly lie at high redshifts z > 1, the most 
straightforward method to increase the correlation amplitude would be to have a strongly 
peaked redshift distribution. Such a peaked distribution, though in principle possible is 
unlikely given the broad shape of the luminosity function, unless special galaxy formation 
models are invoked, such as by Babul and Ferguson (1995). Furthermore, the redshift 
that correspondes to the 15 percentile point in the high redshift tail of a standard galaxy 
distribution is about z = 0.76 at our magnitude limit. Any subset of this tail will not have 
enough galaxies to make up the observed 15%. 

Another explanation is that the very blue galaxies are predominantly a rapidly merging 
population at z > 1. In such a scenario, however, one might expect these galaxies to have 
a correlation function with an excess of close pairs. While several studies support an 
increase of close pairs at the level of (1 + z) 3 (Zepf & Koo 1989; Burkey et al. 1994; 
Carlberg, Pritchet, & Infante 1994; Yee & Ellingson 1995), the fraction of participating 
galaxies remains small, at the level of 10% or less. More importantly, very close pairs do 
not appear to have colors that differ significantly from typical field galaxies. Since we also 
find that the very blue galaxy correlation function has the same smooth power-law slope as 
that of the red galaxies, without a break to higher amplitudes on small scales, we consider 
the merger hypothesis an unattractive explanation. In any case, it is not the existence of 
close pairs that accounts for the enhancements in our signals. 

The other extreme possibility is that the blue galaxies are at low redshift. Indeed, the 
enhanced correlations suggest this to be the case — in the integral that gives the angular 
correlation amplitude, the (dn/dz) 2 term has a very large weight at low redshifts, thus it is 
much easier to account for the excess clustering by modifying the redshift distribution at 
the low z end. To quantify this result, without going into an overly complicated modeling 
scheme, we will try to parameterize the phenomenological redshift distribution of the blue 
galaxies by a simple functional form. This will show that, independent of the particular 
details, the only possible solution is to have most of these galaxies at redshifts z < 0.3. 

For this simple parametrization we make the following assumptions: (1) the number 
counts combined with redshift surveys fix the redshift distribution of the whole sample, 
with 20 < Bj < 23.5, (2) the red galaxies in the sample do not evolve, and their redshift 
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distribution is adequately described by the models, (3) a standard, non-evolving clustering 
amplitude is used throughout. 

For our first experiment, we try to assign the blue galaxies to as low a redshift as 
possible. We adopt the following algorithm for selecting the redshift distribution of the 
blue subset: we subtract the red galaxies from the overall redshift distribution, creating 
a complementary sample. This subset is multiplied by a simple lowpass window, and the 
resulting distribution is assigned to the bluest subset. The overall population fraction 
and amplification of w{9) of this subset are calculated. Satisfactory subsets are required 
to have a fraction greater than 15%, and an amplification greater than 6. Although this 
recipe falls short of providing a full modeling of all the color distributions, number counts, 
and correlations, it is straightforward and incorporates some basic galaxy distribution 
properties. The window to select the low-redshift galaxies is a simple, sharp low-pass 
filter, i.e. all galaxies below that redshift are assigned to have such blue colors. This is 
clearly unrealistic, but this is how we can make the redshifts of the blue galaxies the lowest 
possible. We find that the highest redshift cutoff with such models is at z = 0.28. 

This simple model is clearly in conflict with the redshift surveys. Glazebrook et 
al. (1995) found that at least 4 of the 21 bluest galaxies were at redshifts beyond 0.6. 
Adopting an additional conservative constraint, that at least 25% of this bluest subset 
has to be at z > 0.8, we can do another experiment. We start with the bluest galaxies 
from the models. Their redshift distribution is quite flat. We add an extra low redshift 
component from the rest of the galaxies, with a sharp upper cutoff in their redshifts. The 
highest cutoff satisfying these constraints is at z = 0.2, in which case 60% of the galaxies 
are below z = 0.3. In this scenario all of the remaining galaxies had to be assigned to the 
low redshift tail. The other extreme case in such a picture is to assign extra galaxies at 
very low redshifts, at z < 0.05. In this case we can get away with adding a small fraction 
of all the blue galaxies; however, the small fraction necessary for the amplification exceeds 
the number of available galaxies at that redshift by a factor of 2. 

5. Conclusion: Low- Redshift Faint Galaxies 

We conclude from these exercises that it is not possible to explain the observed excess 
correlations with even moderately large modifications of the class of models we have con- 
sidered, that is, a substantial fraction of the blue galaxies must exist below a redshift of 
z = 0.3. However, such analyses alone cannot discriminate between a population of blue 
galaxies all at redshifts up to 0.3 (e.g. low surface brightness galaxies), versus a strongly 
evolving and thus more peaked population. However, we can certainly rule out the pos- 
sibility that the bulk of these galaxies are of high luminosities with a wide spread over 
high redshifts. From the results of the data and the models, we conservatively predict that 
more than 50% of the currently unidentified population of the faint blue galaxies at our 
magnitude limit will be found at redshifts z < 0.3. 

With evidence that these very blue galaxies are mostly nearby, their luminosities can 
be estimated by applying a K-correction of 0. The galaxies in this survey typically have 
22 < Bj < 23.5, which translates to limits on absolute magnitudes of —18 > Mb,, > —19.5 
at redshifts of 0.3. For lower redshifts, the derived luminosities will, of course, be fainter. 
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Typical (L*) galaxies have Mj ~ —21.0 and thus the very blue galaxies appear to be of 
low luminosities. If their blue color is due to enhanced star formation activity, so that 
their luminosity has increased as well, their progenitors may well be from the faint end of 
the luminosity function. 

Since a rapid luminosity evolution of these blue galaxies, combined with our magnitude 
cutoff at Bj = 23.5, can create a reasonably peaked redshift distribution between 0.2 and 

0. 3, a low redshift population would result if the faint end slope of the luminosity function 
of very blue galaxies were steep, as suggested by the results of Marzke et al. (1994) or 
Metcalfe et al. (1991). Such a hypothesis is not only supported by our work, but also 
by recent redshift surveys. Consider the expected color tracks traversed by galaxies of 
different intrinsic colors as a function of redshift in the Canada-France Redshift Survey 
(see for example Fig 2 of Lilly et al. (1995) which shows (V — I)ab y s redshift to i" < 22). 
The exact distribution of these galaxies involves a complex interplay between the depth of 
the survey, luminosity functions, evolution, cosmology, dust; however, the bluest galaxies 
are found to be at either very low or very high redshifts. Based on this figure of almost 
600 galaxies, we find that the bluest galaxies indeed appear to be concentrated at low 
redshifts. More specifically, reading off the uncrowded bluest 18 points, 11 are below 
z = 0.3, just as we have predicted with our correlation analysis. The strongly evolving 
luminosity function of the bluest half in this data set also suggests that these effects are 
even more pronounced for the bluest 15%. Therefore, we conclude that the faint end 
slope of the local luminosity function for blue galaxies is steeper than that adopted by our 
baseline model and that perhaps the luminous component at high redshifts is less. Thus 
this redshift survey supports our work and shows the effectiveness of our approach. Of the 
10 bluest galaxies with redshifts in the Glazebrook et al. (1995) sample, 4 had redshifts 
below z = 0.3, and 4 were above 0.6. The remaining 11 are a large enough percentage that 
their result may still be consistent with ours. We predict that these remaining objects will 
be found to be low-metallicity systems with redshifts less than 0.3. 

Having identified the blue galaxies as being dominated by a low redshift population, 
we can estimate the expected strength of the cross-correlation between the red and blue 
subsets. From the models - with the assumption that the intrinsic clustering is average - 
we obtain expected amplifications between 0.7 and 1.3, which yields an amplitude range 
of 18 ± 15. Contrary to this expectation is the observed value —18 ± 10 shown in Table 

1, a discrepancy of 2a. However, our null hypothesis of average clustering is likely to be 
incorrect; unfortunately no local information is available on the intrinsic, spatial cross- 
correlations of different galaxy types. This measured amplitude is marginally consistent 
with a null signal. If the numbers are taken literally, we may speculate that the weak anti- 
correlations indicate that only a small fraction of the very blue galaxies can co-exist with 
the red galaxies at similar, moderate redshifts, and for those that do co-exist, they occupy 
separate environments. Since red, early type galaxies populate high density regions, the 
faint blue galaxies should be in surrounding, segregated lower-density regions, which is not 
inconsistent with common beliefs. The primary limitation in this analysis, which prevented 
us from being able to make much stronger statements about the redshift distribution of 
these faint galaxies, is the poor knowledge of the local auto-and cross correlation properties 
of the different galaxy types. 
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We have shown that angular correlations as a function of color can serve as a very 
sensitive statistical tool to study galaxy evolution as well as to provide new insight into the 
nature of galaxies of extreme colors. Since multicolor photometric observations easily reach 
beyond the practical limits of reliable redshift measurements, angular correlation functions, 
combined with color cuts, will play an increasingly important role in understanding faint 
galaxies in the future. As the local galaxy population is better understood, the accuracy 
of this method will increase substantially. 
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Fig. 1) Results for A w as a function of limiting color assuming a power law 

slope of -0.8 are shown, integrated between 9 and 144 arcseconds for both fields. 
The signals are calculated using all galaxies less than or greater than a given 
U-Rp color for the blue and red respectively. Error bars are Poisson in the 

bin counts (although see Section 4.1). The dashed line shows the result for the 
combined data sets with no cut in color. Notice the enhancements over a factor 
of 10 for both the extreme blue and red subsets. The fraction of the entire data 
set used at each limit is shown below each plot in b&d. 




Fig 2) Results for A w with the data sets divided into bands of 

U-R F color of one magnitude width, assuming a power law slope 

of -0.8, and integrated between 9 and 144 arcseconds for both fields. 
These results show that the enhancements in both blue and red 
are robust over 1.5 magnitudes in U-R p The dashed line shows the 

result for the combined data sets with no cuts in color. 
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Fig 3) The angular correlation function W(0) for U-R F >1.5 (Red) and U-R F <-0. 

(Blue) for both data sets broken up into 5 bins of equal log width between 
18 and 289 arcseconds along with a fudicial line for a no-evolution model with 
q o =0.01 and H Q =50 km sec^Mpc" 1 . The results for both fields in 
both colors are consistent with the canonical power law slope of -0.8. 



